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Abstract 

Critical behaviors of quark-hadron phase transition in high-energy heavy-ion colli- 
sions are investigated with the aim of identifying hadronic observables. The surface of 
^ , the plasma cylinder is mapped onto a 2D lattice. The Ising model is used to simulate 

' configurations corresponding to cross-over transitions in accordance to the findings of 

'nJ" ' QCD lattice gauge theory. Hadrons are formed in clusters of all sizes. Various measures 

, are examined to quantify the fluctuations of the cluster sizes and of the voids among the 

I clusters. The canonical power-law behaviors near the critical temperature are found 

■ for appropriately chosen measures. Since the temperature is not directly observable, 

attention is given to the problem of flnding observable measures. It is demonstrated 
that for the measures considered the dependence on the final-state randomization is 
weak. Thus the critical behavior of the measures proposed is likely to survive the 
scattering effect of the hadron gas in the final state. 
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1 Introduction 



> 

X 

■ The subject of this talk can be introduced by commenting on the four conceptual items 
mentioned in the title, but in reverse order. First, high-energy nuclear collisions. The goal 
is to create quark-gluon plasma. So far the conventional signatures are aimed at detecting 
the primordial quark matter. Suppose that is found. Then what? If we can learn from 
the condensed-matter physics, then we know that critical phenomena form a vibrant area 
of research. That leads to the second item in the title: critical behavior. Our concern will 
be the quark-to-hadron phase transition (PT). The hadron-to-quark PT is, in my view, too 
turbulent in heavy-ion collisions to lend itself to systematic studies, despite some recent 
attempt to use percolation as a model to characterize it. The third item in the title is 
"measures." For a thermal system the heat capacity can be a measure that exhibits the 
critical behavior C oc |T — Td ", while for a magnetic system the magnetic susceptibility 
behaves as x |T — Td"'^. What is the measure for quark-hadron PT? Is there a measure, 
call it /i, which behaves as /i ~ |T — Tc|~'^? If so, how do we find it? That brings us to the 
last item: observable. Temperature is not directly observable. (The transverse momentum 
Pt is an unreliable, gross indicator of the average phenomenon.) Moreover, the system is 
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not static and the emitted hadrons may undergo final-state interactions that can smear the 
signature that is to be observed. The aim of this talk is to show how all these issues are 
addressed and to indicate where the answers are to be found. 

The measure that I believe will convey the information about quark-hadron PT is about 
fiuctuations, not only bin-to-bin fiuctuations, but also event-to-event fiuctuations. A way to 
illustrate this belief is to consider the annual precipitation (AP) of two places with widely 
different climatic characteristics. The AP of certain regions in Oregon and Florida are very 
similar, about 40" per year, but their rainfall patterns are nearly opposite: drizzly in Oregon, 
stormy in Florida. The problem is that in determining the PA the daily rainfall is integrated 
over the entire year, and then averaged over all years. That is exactly how the average 
multiplicity of heavy-ion collisions is determined: integrated over the whole hadronization 
time of an event (which can take over 30 fm/c), and then averaged over all events. Clearly, 
to learn about the nature of rainfall or the properties of hadronization, it is necessary to 
study spatial and temporal fluctuations. 

2 The Problem 

If quark-gluon plasma is indeed formed in a heavy-ion collision, it is reasonable to assume that 
the plasma occupies a cylindrical volume, expanding rapidly in the longitudinal direction, 
and more slowly radially. The medium most likely has a temperature profile that is high 
at the center of the transverse plane and decreases with increasing radius due to the radial 
expansion and cooling. When T reaches at the surface, the quark matter undergoes a PT 
into hadrons. Thus the characteristics of critical behavior are to be found on the cyhndrical 
surface at any given time of the evolutionary history of the system. That is like raining 
at a particular instant of observation. If we collect all the hadrons emitted throughout the 
long hadronization processes, the features of critical behavior are likely to overlap and get 
averaged out. It is therefore important to select a narrow At cut in the hadronization time. 
If the emission time and the transverse momentum pt of the hadrons are correlated in a one- 
to-one relationship, then a narrow Apx cut would accomplish what we need. Unfortunately, 
that relationship is not one-to-one. On the contrary, at each emission time the emitted 
particles have a very wide Pt distribution. That wide distribution turns out to provide us 
with a way to proceed. Suppose we make a narrow Apx cut and study the hadronic patterns 
in the rj-cj) 2D space orthogonal to px- Then at any emission time only a small fraction of 
the emitted hadrons would enter the narrow Apr window. This random selection process 
at each time frame thus decorrelates the hadrons collected over the entire hadronization 
process. The measured hadronic patterns in the ri-(f) space is then the superposition of small 
parts of many configurations (which we call configuration mixing), each of which we shall 
stimulate before mixing. The object of our analysis is to find the scaling property of each 
configuration that exhibits the critical behavior, and then to see whether it survives the 
mixing process. 

It is well known that a system at the critical point forms clusters of all sizes that exhibit 
Kadanoff scaling. The details of the dynamical system are not important, since the critical 
behavior depends mainly on the dimension and symmetry of the system. On the basis of that 
universality we shall use the Ising model as a simple device to simulate the configuration 
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of hadronic clusters formed on the surface of the plasma cylinder at each time frame of 
hadronization. There are reasons to believe that the QCD dynamics for the chiral PT 
suggests a cross-over near the second-order PT that belongs to the same universality class as 
the Ising system We shall define hadron formation on the Ising lattice in such a way as to 
be able to simulate a cross-over. From the simulated configurations we can see the hadronic 
clusters randomly scattered over a background of voids. In our first-attempt analysis of the 
patterns formed, we shall ignore the complications that arise from the px distribution and 
consider directly the scaling properties of the clusters and voids in uncorrelated configurations 
that we shall simulate. 

3 Cluster Analysis 

The use of the Ising model on a 2D lattice to simulate second-order PT is standard 0, ^. 
To have a cross-over in quark-hadron PT without the use of an external field in the Ising 
Hamiltonian is not standard [Q. The method uses a cell on the lattice, of dimension e x e 
where e may be taken to be 4, to define hadron density pc at the cth cell 

p, = A02^(0,) (1) 

where 

j&C j£L^ 

The first sum in (0) is over all sites in a cell having site-spins aj; the second sum is over 
all lattice sites, giving the total magnetization m/,. We use the direction of rriL of each 
configuration to serve as the direction of an external field and define the order parameter 0c 
as the cell spin at c relative to that direction. A in (0) is a parameter relating hadron density 
to Ising spins; all measurable quantities should be insensitive to A. Because of the 6'-function 
in (|l|), Pc is either positive where hadrons are emitted, or zero where a void exists. 

There are two quantities we have found that can effectively serve as measures of fluctu- 
ations of hadronic clusters from bin to bin and from configuration to configuration. One is 
J(M, T) defined 

and the other is K(M,T) whose definition is omitted here, but can be found in Ref. 0. In 
(H) (■ ■ ■) refers to averaging over all bins [M of them, M = {L/5)^] and over all configurations. 
We discuss only the properties of J{M, T) here; those of K{M, T) are similar. 
What we have found about J{M,T) is that it is factorizable p, i.e., 

J{M,T) = a{T)Jc{M) , (4) 

where Jc{M) = J{M,Tc). Thus, by definition, a (Tc) = 1. Yet, when T is away from T^, but 
not too far in T < T^, a;(T) satisfies the power-law behavior 

a(r) oc (Tc - r)-^ C = l-88 . (5) 
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This may be regarded as the critical behavior that we have been searching for. However, T 
is not directly measurable. We can use p instead, where p is the cell density averaged over 
all cells in one configuration. Then we find 

J{Mrp) = a{p)UM) (6) 

with 

a(p)(x(p-po)"S C = 0.97 , (7) 

where po is the average density where J{M,po) = Jo{M) is the highest. The behavior in (0) 
is measurable and can be used to signal the fiuctuation property of a quark-hadron PT. 



4 Void Analysis 

In addition to finding the critical behavior of the hadronic clusters, we can also analyze the 
scaling properties of the voids where no hadrons appear @]. First, it is necessary to define 
precisely what a void is, and how to avoid a single hadron from changing its quantification. 
To that end we define a bin to be empty if the average density pk of the kih bin is less than 
Po, where po is a parameter that can be varied by the analyst of the data. We do not set po 
to zero in order to eliminate the small fiuctuations due to irregular occurrences of hadrons 
in or at the edges of a bin. A void is then the area of contiguous empty bins, and its area is 
therefore 

Vn = J2HP0-Pk) , (8) 

k 

where the sum is over all connected empty-bins of the nth void. Vn is in the unit of number 
of bins. A quantity less dependent on the size of the bins is the fraction Xn = Vn/M. For 
every configuration, there exists a set of void fractions {xi,X2,X3, . . .Xm} where m is the 
total number of voids. The comparison of void patterns of different configurations would be 
difficult if they are characterized by those sets. To ease that comparison, we define the G 
moments 

W 

where {x'^) = (Z]n=i ^1) /"^' average performed for each configuration. Thus Gg is a 
number (for each q) that characterizes the void pattern of a configuration. At the critical 
point this Gq fiuctuates widely from configuration to configuration. It depends on the bin 
size. 

At T in the vicinity of we find that (Gg), averaged over all configurations, satisfies 
scaling law for a variety of po values, i.e., 

(Gg) oc M> . (10) 

That means voids of all sizes occur. Furthermore, the scaling exponents depend on q linearly 

7<7 = Co + cq (11) 

with c = 0.8 at Tc and po = 20. This value of po is less than 8% of the maximum p. The value 
of c depends on T and po mildly, and can be used to quantify the nature of void patterns. 



G, = f^ , (9) 



4 



Since Gq fluctuates very widely, i.e., the distribution P{Gq) is very broad, the moments 
of Gq can reveal that fluctuation. We have studied the derivative of the first moment, which 
is 

Sq = {GqinGq) (12) 

and found that it also satisfies a scaling law 

Sq oc (13) 

whose exponent depends on q linearly 

(Tg = So + sq (14) 

with s = 0.76 at T = Tc and po = 20. The value of s also depends on T and po niildly, and 
therefore can further be used to quantify the nature of the fluctuation of void patterns from 
configuration to configuration. 

The dependences of c and s on po and T, though mild, can be used to diagnose the 
temperature range of quark-hadron PT. For po > 20 the values of c and s determined in our 
analysis vary in a range as much as 20% for T varying in the range of ±2% around Tc, but 
they converge to the values quoted as po 20. If c and s are not unique for a range of T 
around T^, it means that jq and are not fixed (for fixed q) if PT takes place over a range 
of T, not exactly at Tc. Then there would be no scaling behavior. Thus, if the experimental 
data reveal the loss of scaling as po is increased, it implies that the PT takes place over a 
range of T. If, however, scaling persists, then PT occurs only at one T, presumably Tc. 



5 Complications in Hadronization 

There are two aspects about hadronization in heavy-ion collisions that one should be con- 
cerned about for the type of measures we are proposing. One is the effect due to final-state 
interaction as the hadrons traverse the hadron gas. The other is the effect of configuration 
mixing. The former is studied by requiring the produced hadron in each configuration to 
take u steps of random walk, since scattering by a hadron gas implies the random defiections 
of the initial pr that result in random wandering in the r/-^ plane. We have studied the 
effect of such random walks (for u up to 6) on J(M, T) and found that the factorizability 
(|) persists and that a{T) retains the essential character of (|^) 0. 

The other complication of configuration mixing due to the non-correspondence between 
Apx and At is more difficult to treat. We have considered only the simple case of no corre- 
lation between successive time frames. We mix four independently simulated configurations 
by taking one quadrant from each and then apply our void analysis to the mixed configu- 
rations 0]. We find that the effects on 7^ and aq are essentially negligible. This problem 
will be followed by a more realistic treatment of the configuration mixing where the wide pt 
distribution of the emitted particles at any given time will be taken into account together 
with the hadron-void correlation in successive time frames. 
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6 Multifragmants 



Finally, let me make some remarks about "Common Aspects" that is a part of the title 
of this session. Despite the valiant attempts made by the organizers of this conference to 
promote interaction between the two subfields encompassed here, cross-fertilization has not 
been a vibrant phenomenon that has taken place. However, it seems that the ideas about 
the void analysis discussed here can be translated into the study of multifragmentation 
in the intermediate mass range. Let me therefore venture a suggestion on an analysis of 
the fragments with the hope that it may yield a feature of the experimental results more 
challenging to the model builders than what they have hitherto dealt with. 

Suppose that for each event an experiment can measure a certain number of fragments 
whose total mass may be less than the total mass of the nuclei in collision and may fluctuate 
from event to event. Nevertheless, for all the fragments Mj that have been detected, one 
can calculate the fraction Xi = Mi/ J2i If it is the charges of the fragments that are 
measured, then let Xi = Zi/J^i^i- Each event is then characterized by a set of fractions 
{xi,X2,X3, . . .Xn}- The idea is to go beyond studying the averages and examine the nature 
of the fluctuation of the fragments detected. Let us define 

1 " 

H,{n) = -Y.xl (15) 

where n is the total number of fragments measure in an event. Clearly, Hq = 1 and Hi = 1/n. 
At higher q, Hq{n) are progressively smaller, but are increasingly more dominated by the 
large Xi components. 

If P (Hg) is the distribution of Hg after all events are sampled, then the average is 

W = ^ E = J dHg Hg P{Hg) (16) 

where e labels an event and J\f is the total number of events. (Hg) provides some average 
information about the multifragments, but to study the fluctuation from event to event one 
should calculate the normalized moments of moments 

Cp,, = (H'g) / {Hgf (17) 

This Cp^g may contain too much information. The derivative at p = 1 may already be 
sufficient to provide some interesting insight into the erratic nature of the fluctuation 



' dp 



W {Hg) 



the study of which may be called erraticity 0, for convenience. The experimental determi- 
nation of Eg as a function of q may pose a serious challenge to the theoretical models, some 
of which may not contain enough dynamical fluctuations to generate the erraticity revealed 
by the data. 
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